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Abstract 


The question whether one may identify a neutron star by observing 
hydrogen and helium lines in its spectra is examined. An analysis is presented 
of the burning and diffusion of hydrogen and helium in a neutron star atmosphere. 
An initial two-component atmosphere of hydrogen and carbon and of helium and 
carbon is examined utilizing various photospheric temperatures. The hydrogen 
bums through the proton-proton reaction and the CN cycle and the helium by 
the triple alpha process and alpha capture on carbon. It is found that only 
a negligible proportion of hydrogen or helium can exist after a few years in 
the neutron star atmosphere. Therefore spectra which include hydrogen and 
helium lines cannot be attributed to neutron stars. 
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I. INTRODUCTION 

There is at present a great deal of uncertainty as to the chemical 
composition of a neutron star atmosphere. This is regrettable since it 
disallows one to associate with any degree of certainty an optical 
source with a neutron star. Work has been performed on the statistical 
equilibrium composition of neutron star atmospheres (Tsuruta and Cameron 
1965). However, because of the extremely small scale heights in the 
neutron star envelope, diffusion effects may drastically change the 
composition. The density changes by about eight orders of magnitude 
and the temperature by about two orders of magnitude in a depth of 
approximately one to twenty meters depending upon the choice of model 
(Morton 1964, Tsuruta and Cameron 1966). Chiu and Salpeter (1964) have 
estimated that hydrogen and helium will not be present in the envelope 
of a neutron star whose age is greater than a few years. The hydrogen 
and helium will have diffused down to the interior regions and have been 
completely burned due to the high temperatures prevailing there. It is 
the purpose of this work to examine this idea in detail: to analyse the 
degree of burning and diffusion of hydrogen and helium in a neutron 
star with the end in mind of noting whether any optical sources contain- 
ing hydrogen or helium lines may be associated with neutron stars . 

II. APPROACH TO THE PROBLEM 

A typical neutron star model of 1.01^ and 10 km. radius has been 
chosen for this problem Initially, the atmosphere of the star is 
composed entirely of carbon with a layer of hydrogen or helium super- 
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imposed on the carbon for the two cases. This initial configuration is 
chosen in order to maximize the time required to burn any fraction of 
the original mass of hydrogen or helium. The hydrogen burns through a 
combination of the proton-proton reaction and the CN cycle; the helium 
reacts through the triple alpha process and alpha capture on carbon. 

Carbon burning is neglected in this study. Although one may expect a 
neutron star atmosphere to contain nuclei up to the iron peak (Tsuruta 
and Cameron 1965, Cameron 1965), the proton and alpha capture rates will 
be lowest for carbon at the relevant temperatures. At a temperature of 
1.0x10^ °K, < ctv > for the CN cycle is 1.88x10 23 cm 3 /sec (Reeves 1965) 

whereas < ctv > for the reactions: 23 Na(p,Y) 2l+ Mg, 31 P(p>y) 32 S and 

07 28 -21 2 -22 2 

Al(p,y) Si is given by 1.02x10 cnr/sec, 2.93x10 cnr/sec and 

4-24x10 cnr/sec respectively (Truran, Hansen, Cameron and Gilbert 

1966). At this temperature, < Ctv > for the reaction '*' 2 C(a,Y)‘*'^0 is 

3.66x10 3 ^ cm 3 /sec compared to 1.51xl0 -2 ^ cm 3 /sec, 1.30xl0~ 2 ^ cm 3 /sec 

and 2.28x10 2 ^ cm 3 /sec for the reactions '*'^0(Q;,Y) 2 ^Mg, and 2 Svig(a,Y) 2 ^Si 

respectively (Truran et al 1965). Therefore the presence of advanced 

products of evolution in the neutron star atmosphere would speed up the 

burning of hydrogen and helium. 

Since the particle diffusion velocity is proportional to the gradient 
of the number density, an intermediate, continuously varying strata of 
carbon plus hydrogen or helium is inserted between the purely carbon 
atmosphere and the overlying layer of pure hydrogen or helium. The 
diffusion takes place under the influence of external forces. These 
forces include gravity and the electric field due to the charge separation 
of ions and electrons in the atmosphere. 
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The time -dependent coupled set of equations of stellar structure, 
nuclear reactions and diffusion are solved numerically by computer. 

The entire atmosphere is divided into zones and the equations solved 
for each zone at each epoch. The resulting variables are then used to 
solve the equations for each zone at the next epoch. This procedure is 
continued until at least .95 of the original mass of the hydrogen or 
helium layer has been burned. 

III. EQUATIONS OF STELLAR STRUCTURE 

The non-equilibrium stellar structure equations for a spherically 
symmetric diffusing system are (Chapman and Cowling 1952, Chandrasekhar 

1939) 
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P = P(p,T) 


(5) 


where P, M, p, T, L and K are the pressure, mass, density^ temperature, 
luminosity and opacity respectively, p^ and c^ are the density and diffusion 
velocity of the i^* 1 component. The radiative energy transfer equation ( 3 ) 
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is utilized since the neutron star is stable against convection (Morton 
1964, Tsuruta and Cameron 1966). 

Because of the special nature of neutron star atmospheres, certain 

simplifying approximations to the equations of stellar structure may be 

made. Previous studies of equilibrium neutron star atmospheres have 

shown that the envelope is only about one to ten meters thick (Tsuruta 

1964, Tsuruta and Cameron 1966). The mass in equation (l) may then be 

set equal to the total mass of the star, thereby decoupling the mass 

equation (4). Further, energy generation in such a thin atmosphere 

is very small compared with interior processes so that the luminosity 

in equation (3) may be taken as a constant. For a neutron star of one 

GM 10 

solar mass and a radius, R , of 10 km. , g = =-2 w 10 • Therefore 

the last two terms in equation (l) are extremely small compared with 
the first two so that a "quasi-static" equation may be used Since the 
depth of the atmosphere constitutes less than one per cent of the radius 
of the star, the radial distance in equation (l) is taken as constant 
and set equal to the radius of the star. Equations (1,3) then reduce 
to 


dP 

dr 


- Pg 


(6) 


dT _ yp 

dr = ' ^ (7) 


where B is a constant equal to 


3L 

l6rrac 


5 


The equation of state, represented as a function of density and 
temperature in equation (4), is that for a perfect gas of electrons and 
ions plus radiation pressure. 


P = P. + P 
i e 
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(8) 


where P^ and P g are the pressures of the ions and electrons respectively 
and a is the Stefan-Boltzmann constant. The general expression for the 
pressure of a fermion of spin 1/2 is given by (Landau and Lifshitz 1958) 


P = 


8tt 

3h : 


; (E-U)/kT + 1 dp 


SE , 
^ d P 


(9) 


where U is the chemical potential of the particle. The contribution to 
pressure by the ions is only important in regions in which degeneracy 
effects are unimportant . Therefore the ion pressure may be represented 
throughout the atmosphere by the non-degenerate perfect gas equation of 
state 


% sl i H = f ( 10 > 

where H is the mass of a proton and = p/(nJi). 

It is necessary to evaluate the integral in equation ( 9 ) for the 
electrons, however, at each density and temperature in the atmosphere. 
This is done most easily by tabulating the equation of state for the 
electrons in terms of a degeneracy parameter \ (Handbuch der Physik 


1958) 
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( 12 ) 


where the functions anc ^ F^y^X) are defined by 


n oo 


F s (X) = 


t s dt 

t-X 

e +1 


(13) 


and (Jt g = p/(n g H). In the non-degenerate regions of the envelope, the 
ratio(2/3)F^/2(^)/F-|y2^) approaches unity and the equation of state 
(ll) reduces to the non-degenerate perfect gas law. As the density 
increases for a given temperature, the parameter X increases. When X 
exceeds 20 , the electron gas may be considered to be completely 
degenerate (Handbuch der Physik 1958). At this point, the Fermi energy 
of the electron is sufficiently high compared with kT such that a zero 
temperature Fermi gas well approximates the statistical distribution 
of electrons. The equation of state is then given by (Chandrasekhar 


1939) 


where 


li c; 
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and 
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IV. PARTICLE DIFFUSION 


A two particle diffusion approximation will be utilized in this 
calculation to find the diffusion velocity of hydrogen and helium in 
the neutron star atmosphere. The relative diffusion velocity of such 
a two particle gas mixture is given by (Chapman and Cowling 1952) 



— — p 

n L n 2 v 12 




(18) 


where n = n^ + n 2 , D^ 2 is the coefficient of diffusion, is the thermal- 

— » 

diffusion ratio and d 12 is defined for a gas with zero mass flux to be 

^12 = iiKT L K V ( n l T ) " P l^l/ m l] 


Fj_is the resultant external force acting on species 1. For the problem 
under consideration, this includes the gravitational force and electric 
field due to the charge separation of ions and electrons in the gravita- 
tional field. 

Both the coefficient of diffusion and the thermal-diffusion ratio 

depend on the form of the particle interactions which in turn are a 

function of the degree of ionization in the envelope. The temperature 

in the atmosphere of a neutron star depends upon the model chosen and 

on the time which has elapsed since formation of the star. Tsuruta and 

Cameron (1966) have shown that for neutron stars of about one solar 

3 5 

mass, it takes approximately 10 to 10 years for the photospheric 
temperature to cool from 10^ °K to 10^ °K. During this time, the neutron 
star may be observed as a point source black-body emitter with a 
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maximum in the x-ray region. Since the significant temperature range of 
the envelope is at least 10^ °K, complete ionization of the carbon- 
hydrogen-helium atmosphere will be assumed in this calculation. Thus 
the particle interaction is that of Coulomb repulsion between completely 
ionized particles and the coefficient of diffusion is given by (Chapman 
and Cowling 1952) 


D 




ln[l + ctn(9 0 / 2 )J 


( 20 ) 


and Zg are the charges of the two ions 
and 0 q is the minimum scattering angle. Chapman and Cowling use the 
average for ctn(9 Q /2) given by 


where iy = nyy/tny + iy), Z ± 


ctn(e o /2) 



(20a) 


where d is the mean interparticle-spacing of ions taken in this work to 
be [3/(^ n )] 1//3 - It is quite possible that this underestimates the 
number of small angle collisions thereby giving too large a diffusion 
coefficient and too large a relative diffusion velocity. It is the wish 
of the author to maximize the time required for diffusion and burning 
whenever a choice of variables may be made. Therefore, in accordance 
with the techniques of Mestel (1949 ) 5 Lee (1950) and Spitzer (1956), 
the minimum scattering angle 9 q is defined as 

e _ \ a = h/CBkTp/n] 1 / 2 ^ 

0 ^ D [kl/(4TTe 2 )] 1 / 2 /(Z 1 2 n 1 + Z^n 2 ) 1 / 2 (20b) 

where is the de Broglie wave length and X^ is the Debye 
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screening radius for the two ion mixture. The method employed in this 
analysis is to calculate ctn(9 o /2) through both equations (20a) and 
(20b) and choose the largest value in order to minimize The 

contribution of Vt to the diffusion equation (18) is small for a 
neutron star and will be neglected (Chapman and Cowling, Tsuruta and 
Cameron 1966). 

Since the time scales for diffusion are very slow compared to the 
mean collision time of the gas mixture, the distribution of particles 
will differ only slightly from statistical equilibrium. Therefore the 
force will be taken as that for the equilibrium configuration and is 
given by 

F 1 = - m^gr + Z^eE (21) 

2 

where g = MG/r^ as in equation (6) . The gravitational field produces 

a separation of electrons and ions which produces the electric field 
— > 

E in equation (21). This field, which acts radially outwards and tends 
to reduce the effect of gravity on the ions, is given by (Spitzer 1956) 




The electric field reduces the magnitude of the gravitational force by 
a factor of l/(Z-^ + l). 
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IV. NUMBER DENSITY NETWORK 


The number density of a particular species can change in two ways: 

by diffusion and by nuclear burning. This may be represented for the 

.th . 

1 species by 


dn. i -v , 

rr = -s h l 


2 \ 

n. c. J 

i i/ 


(24) 


where Q(n^) is the nuclear burning and production rate of species i. 

Hydrogen burns to produce helium through the proton-proton reaction 

and the CN cycle and helium burns to produce carbon by the triple alpha 

12 l6 

process and through the reaction C(a,Y) 0. At any given time and 
position in the atmosphere, the number density network for the three 
nuclei: hydrogen, helium and carbon is given by 



(25) 


dn He 
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r 

The reaction probabilities; 
(Reeves 1965 )* 


( 26 ) 

) + ! V 3 + f Vc (27) 
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S, Y> & and e (in cm 0 / sec) are given by 
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P = ( 1 . 0+2 . 5P 1 / 2 /t 6 3 / 2 ) (1 . 0+1 . 23xl0" 2 Tg i//3 +7 • 8xlO _3 Tg 2 / 3 +6 . 73x10"^ ) 

X (3.3x10" 13 h) Exp ( -33 - 8l/T 6 1 / 3 )/t 6 2 / 3 (28) 

Y = ( 1. 0+1. 75P 1 ^ 2 /T 6 3 / 2 )( 1.0+2. 74x10" 3 T 6 1 / 3 - 3-74 x10' 3 T 6 2 / 3 -7.17x10 _5 T 6 ) 

X (4.79x 10 9 H) Exp (-152.3l/T 6 1//3 )/T 6 2 / 3 (29) 

6 = Exp (2.4 p 1 / 2 /t 6 3 / 2 ) H 2 [9.47x10" 6 • 10 " 

/ 119.0 % . 

+ 5.0x10 5 • 10 K Tq ’ ]/ Tg° ( 30 ) 

e = Exp(2.4p 1 / 2 /T 6 3//2 ) • H* {[l.84xl0 10 • io"^ 3 °'° 5 / T 8 ^ j 

Tg 2 (l-4xl0" 3 Tg + 0.2Tg" 2 / 3 ) 2 ] + ^4.3 x 10 3 .10" 122//T 8 + 2 •35xlO i+ xlO" 135 / T 8j / /Tg 

2/l 1/3 ^ 3 " L ) 

6 = 0.06Tg /5 (1 + 0.07Tg 7 J ) (31a) 

where Tg = T/l.QxlO^ °K and Tg = T/l.GxlO^ °K. The helium which is 

burned through the reaction ^" 2 C(ck,y)^ 0 is recycled into carbon. This 

is done to simplify the atmosphere to a two-component one for helium 

burning and more important to maximize the burning time. Reeves ( 1965 ) 

has shewn that for temperatures greater than 3-4x10^ °K, the reaction 

’'"^ > 0(Q',Y) 2 ^Ne will always be faster than the reaction ''" 2 C(af , y)^ 0. Since 

all significant helium burning through alpha capture will take place 

above this temperature, the procedure utilized in this work will slow 

down the total time required to burn the helium. 


V. COMPUTATIONAL PROCEDURE 


The solution of the diffusion and burning of hydrogen and helium 
in an atmosphere of hydrogen, helium and carbon is a time-dependent 
problem which requires coupling the equations of stellar structure, 
equations (6, 7, 8), with the diffusion equations (l8, 19) and the number 
density network, equations (25, 26, 27)- This is accomplished with 
the use of a computer by differencing all equations, choosing boundary 
conditions at the photosphere of the neutron star and then at each 
time step solving for the variables downwards from the photosphere until 
the isothermal internal core temperature is reached. The process is 
then repeated for the next time step utilizing the variables computed 
from the preceding epoch. This procedure is continued until the mass of the 
hydrogen or helium is reduced to a specified fraction of its original value. 

The boundary values of density and pressure at the photosphere are 
found in the following way. For a given radius and luminosity of the 
star, the effective temperature is defined by 


T = 
E 



(32) 


where a is the Stefan's radiation constant. Coupling the definition of the 
optical depth T for the medium: = - ?6(p ,T)p (33) 

with equation (5) yields 
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Integrating from t = “ to the photosphere one finds 


zdT 

fp,T) o 


The Eddington approximation gives T =2/3 and P = aT 4 /6 . For a 

Tin ' O E 


constant g and assuming is independent of height, equation (35) yields 


*> h 3*TP ph ^ 


|t 4 

Z E 


The equation of state gives a second relation for the photospheric 


pressure. 


1 4 P r)h ^ 

p = i aT 4 +,-£2 

ph 3 E (| T + ^ i )H 


Equations (36, 37) may be iterated to find P , and p An initial value 

ph ph 

is chosen for the opacity giving a value of P^ from equation (36). This 

result is substituted into equation (37) to find a p^. This value of 

the density is used to compute a new K which in turn is used to find a 

new P ,. The procedure is continued until the values of P , and p , 
ph ph ph 

converge to within the desired accuracy, in this case one part in 10 4 . 

The opacity is found from a table prepared from the code developed by 
A.N. Cox and his associates of Los Alamos Scientific Laboratory. The 
code includes contributions from bound-bound, bound-free, free-free 
absorption, electron scattering, negative ion absorption and electron 


conduction. 
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Equilibrium studies of neutron star atmospheres have shown that the 
density changes by about 8 orders of magnitude and the temperature by 
about 2 orders of magnitude from those values at the photosphere in a 
distance of less than 20 meters (Tsuruta and Cameron 1966). Because of 
these large changes in the variables with distance, the stellar structure 
equations are written in logarithmic form with pressure as the dependent 
variable in order to reduce inaccuracies for a given zone size. Equations 
(6, 7 ) transform to 

f-§-§ = - | exp (in P - In r - In p) (38) 

^ ^2- £ = exp(ln P + In K - 2 In r - 4 In T) 

d In P 16 tt acg 

(39) 

All equations in this study are differenced implicitly. This is done 
to insure maximum stability in the coupled equations of stellar structure, 
nuclear reactions and diffusion. Severe problems of stability may arise 
should explicit or symmetrical differencing be used with a poor choice 
of time steps. Richtmeyer (1957) has shown that the linear second order 
diffusion equation is unconditionally stable for all time increments 
when differenced implicitly. However, the use of such an approach creates 
difficulties of solution in that the variables to be solved for at the 
updated time appear in all derivatives and in all functions . Thus , no 
simple recursion relations exist as they do for explicit differencing 
techniques between the variables at time t + At and the variables 
at time t. The implicitly differenced equations in this study are 
non-linear and must be solved by some iterative method. The program 
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developed first chooses the proper equation of state which is dependent 
upon the value of the degeneracy parameter X, defined by equation ( 13 ). 
The choice of one of equations (10, 11, or 14) yields a relation between 
the density, temperature and pressure. The substitution of this into 
equations (39> 40) gives two non-linear equations in two unknowns, In r 
and In T at the updated time. For a given In P, these may be written 
as 

p -t + At 

F |_( In T) J+1 , (in I)j , (in r) J+1 , (In r)j j - 0 

(40) 

G [_(ln T) J+1 , (In T)j , (in r) J+1 , (in r)j J =0 

(41) 


where the subscripts represent the zone at which the variables are to be 
found and the superscript indicates that these variables are evaluated 
at the updated time, t + At. The solution of these equations is begun 
with J=1 corresponding to the photosphere. The boundary values are 
therefore assigned to (in T)^ and (in r)^. Subsequent values of the 
variables are found from the values at the preceding zone by utilizing 
the Newton-Raphson method to solve equations (40, 4l). That is for any 
zone, these equations are linearized to yield 


F + c-r: — 6 ( In T ) , + v / 6 ( In r ) T , 

d (In T) J+] _ v 'J+l d (In r) J+] _ v 'J+l 


= 0 
(42) 


G + 5-xrrlr^ 6(1b t 4+i + 6(ln r 4 + i - 0 


'J+l 


( 43 ) 
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6(ln T) t , and 6(ln r) are found for an initial "guess" of (in T) , 

J *X i-L J » -L 

and (in r) T , . New values of the variables are then computed from 
J » X 

(in T) J+1 = (in T) J+1 + 6(in T) J+1 (44) 

(ln r) J+l * ( iri r ) J+ l * 6 < ln r >j+i Ct5) 

These values are then substituted into equations (42, 43) and the process 

is repeated until values of (in T) t ,, and (in r) T ., are found to within 

d+X J+X 

the desired degree of convergence. This process is repeated for 40 
zones. The evaluated temperature and density profile of the atmosphere 
is then substituted into the diffusion and reaction network equations 
to yield new mass fractions at the updated time. For a given initial 
layer of either hydrogen or helium, this process is repeated until at least 
.95 of the mass of the original layer is burned. 

VI. RESULTS AND CONCLUSION 

Figures 1 and 2 illustrate the initial density and temperature profile 
as a function of distance from the photosphere (R - R) of five models 

ili 

employed in this work: two of hydrogen and two of helium, with photospheric 
temperatures of 5 • 0x10^ °K and 1.0x10^ °K and one of hydrogen with a 
photospheric temperature of 2.0x10^ °K. The layers are purely of hydrogen 
or helium for 50 cm. and then vary sinusoidally to zero for the next 100 cm. 
The initial mass fractions for the two component atmospheres are given 
by 


17 


x = 1.0, z = 0.0 R -R ^ A-B (46a) 

hi 

x = sin^ ^ ? z = 1.0 - x A-B ^ Rg-R ^ A (46b) 

x = 0.0, z = 1.0 A ^ r e -R (46c) 

where A and B are 150 cm. and 100 cm. and x and z represent the mass 

fractions of hydrogen (or helium) and carbon respectively. This layer 

will be referred to as the 150 cm. layer. Later in this section a 35 cm. 
layer will be discussed which is described by equations (46a , b, c) with 
A = 35 cm. and B = 10 cm. One notes in Figures 1 and 2 the typical 
behavior of the temperature and density of the neutron star atmosphere: 
the temperature varying appreciably in a distance of less than 10 meters 
and the density changing by four to five orders of magnitude in a distance 
of about 3 meters . 

In Figure 3? the burning and diffusion of hydrogen layers with time 

is compared for the three different photospheric temperatures of 

5.0x10^ °K, 1.0x10^ °K and 2.0x10^ °K for the same thickness (150 cm.) 

layer. The ratio of the integrated mass of the hydrogen layer over the 

entire atmosphere at any time to its integrated mass initially is plotted 

with time. One half of the original masses of 3.28x10^ gm., 4.53x10^ gm. 

and 2.45x10 J gm. remains at 9-6x10 sec., 4.8x10 sec. and 2.5x10 sec. 

for the 5.0x10^ °K, 1.0x10^ °K and 2.0x10^ °K photospheric temperatures 

respectively. The increase in the rate of burning with temperature is 

5/2 

due to two effects: the diffusion velocity increases as the T as noted 
in equation (20) and the nuclear burning rate increases exponentially 
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with temperature. For all three cases 0.95 of the initial hydrogen layers 

5 

are burned in less than 3-5xKr sec. Little difference in the rates of 
hydrogen burning are noted when two different thicknesses of hydrogen 
layers are compared in Figure 4. The 150 cm. layer is compared with the 
35 cm. layer for the same photospheric temperature of 1.0x10^ °K. The 

14 

35 cm. layer corresponds to an initial integrated mass of 1.09x10 gm. 

The two layers of hydrogen are 95% burned in less than 2 *0xKr sec. 

A time study of the mass fraction of hydrogen as a function of 

distance from the photosphere is shown in Figure 5- The decrease in 

mass fraction and the diffusion of hydrogen to points deeper within the 

neutron star is illustrated for the case of the 150 cm. initial layer 

of hydrogen at a photospheric temperature of 1.0x10^ °K. The four times 

4 5 1 5 

of 0.0 sec., 4.0x10 sec., 1.0x10 sec. and 3*4x10 sec. correspond to 

integrated mass ratios of 1.0, 0.518, 0.120 and 0.023. 

The burning and diffusion of the 150 cm. layer of helium is compared 

for the two photospheric temperatures of 1.0x10^ °K and 5*0x10^ °K in 

Figure 6; in Figure 7, the 150 cm. layer of helium is compared to the 

35 cm. layer of helium at the same photospheric temperature of 1.0x10^ °K. 

The thickness of the layer is apparently unimportant as the two layers, 

15 17 

initially of 1.38x10 gm. and 1.55x10 gm. , are seen to burn at almost 

g 

exactly the same rate: 95% burned in less than 8.0x10 sec. Temperature 

difference changes the rate of burning by a factor of less than three. 

17 

The 150 cm. layer, initially 9-60x10 gm. , with a photospheric temperature 
of 5*0x10^ °K is 95% burned in approximately 1.8x10^ sec. 

These results all indicate that 0.95 of the original masses of hydrogen 
and helium are burned in less than one year for photospheric temperatures 
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6 o 

greater than or equal to 5*0x10 K. This is in general agreement -with 

Chiu and Salpeter (1964) who estimated that any hydrogen or helium present 

in a neutron star atmosphere would diffuse inwards and be burned in a 

few years. Since the cooling time of the photosphere of a neutron star 

from T^, = 1.0x10 ^ °K to T^ = 5*0x10^ °K is of the order of 100 years 

(Tsuruta and Cameron 1966), it may be concluded that only a negligible 

proportion of either hydrogen or helium may exist in the atmosphere of 

a neutron star after a few years. Indeed the cooling times presented by 

Tsuruta and Cameron (1966) indicate that a medium mass neutron star, 

initially at an internal temperature of about 10 10 °K will cool to a 

photospheric temperature of 2.0xic/ °K in 10^ - 10^ sec. and to l.GxlO' °K 

in 10 - 10 sec. Since the burning time of hydrogen at T = 2.0x10' °K 

w E 

4 7 Q 

is about 5x10 sec. and the burning time of helium at T„ = 1.0x10 K 

E 

is about 5x10 sec., these photospheric temperatures would appear to 
be the appropriate ones to consider. Further, the analysis presented 
herein probably overestimates the burning time due to the choice of the 
slowest diffusion coefficient and the neglect of proton and alpha capture 
on nuclei heavier than carbon. Hence one would not expect to see any 
evidence of hydrogen or helium lines associated with a neutron star<> 

The author wishes to express his appreciation to A.G.W. Cameron 
and R. White for their many helpful discussions and to E. Tech for his 
assistance with the computer programming. 
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FIGURE CAPTIONS 

Figure 1. Initial density profiles as a function of distance from the 

photosphere, R^,. Five models of 150 cm. initial layers of hydrogen 

or helium are plotted for various photospheric temperatures, T , 

E 

for a one solar mass, 10 fin. neutron star. 

Figure 2. Initial temperature profiles as a function of distance from 
the photosphere, R . Five models of 150 cm. initial layers of 
hydrogen or helium are plotted for various photospheric temperatures, 
T E , for a one solar mass, 10 fin. neutron star. 

Figure 3- Hydrogen burning rate for three photospheric models with the 
same initial depth of layer. The ratio of the integrated mass of 
hydrogen in the neutron star atmosphere over the mass of the initial 
150 cm. layer of hydrogen is plotted as a function of time. 

Figure 4 . Comparison of hydrogen burning rate of two different initial 
thicknesses of layer with the same photospheric temperature, 

rj O 

Tg = 1.0x10 K. The ratio of the integrated mass of hydrogen in 
the neutron star atmosphere over the mass of the initial layer of 
hydrogen is plotted as a function of time. 

Figure 5 . Time history of 150 cm. initial layer of hydrogen in neutron 

rj 

star atmosphere with a photospheric temperature, Tg = 1.0x10 K. 

The mass ratio of hydrogen in the atmosphere is plotted as a function 
of distance from the photosphere, R^,, for four times. 



22 


Figure 6. Helium burning rate for two photospheric models with the same 
initial depth of layer. The ratio of the integrated mass of helium 
in the neutron star atmosphere over the mass of the initial 150 cm. 
layer of helium is plotted as a function of time. 

Figure 7 - Comparison of helium burning rate of two different initial 
thicknesses of layer with the same photospheric temperature , 

rj q 

Tg = 1.0x10 K. The ratio of the integrated mass of helium in the 

neutron star atmosphere over the mass of the initial layer of 
helium is plotted as a function of time. 
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